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A definition of topological phases of density matrices is presented. The topological invariants in 
case of both noninteracting and interacting systems are extended to nonzero temperatures. Influence 
of electron interactions on topological insulators at finite temperatures is investigated. A correlated 
topological insulator is described by the Kane-Mele model, which is extended by the interaction 
term of the Falicov-Kimball type. Within the Hartree-Fock and the Hubbard I approximations 
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phase transition into the topologically trivial Mott insulator. 
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I. INTRODUCTION 

Recently discovered topological insulators (TIs) have 
attracted a lot of interest in condensed matter physics. 1-6 
This novel electronic phase in solids has a band gap in 
the bulk but it can conduct an electric current via gap¬ 
less edge states which are robust against scattering on 
impurities or other weak perturbations of ideal systems.' 
Formation of the metallic edge states is related to the 
nontrivial topology of the ground state which originates 
from a spin-orbit coupling. 8 The TIs are characterized 
by topological invariants which cannot be continuously 
changed unless the single-particle gap is closed.' The idea 
of the TIs has been extended into symmetry protected 
topological phases which possess nontrivial topological 
properties as long as the symmetries are present. 8 The 
existence of the topological gapless edge states was exper¬ 
imentally confirmed by transport measurements 9-12 and 
an angle-resolved photo-emission spectroscopy 10,13-15 in 
many materials, such as, e.g., B^Sez, Pb\- x Sn x Se , or 
HgTe/CdTe quantum wells. 

The existence of the TIs and most of their properties 
can be understood and described within the noninter¬ 
acting band theory. Recently, an influence of electron- 
electron interactions on the topological insulating states 
has attracted a growing attention. 16 Theoretical studies 
are focused mainly on two effects of the strong corre¬ 
lations in topologically nontrivial conditions. The first 
effect is when a non-on-site interaction yields new topo¬ 
logical phases. 17 The second problem is how strong cor¬ 
relations change properties of the band topological insu¬ 
lators. Namely, interaction can lead to the creation of an 
ordered phase which through renormalization of the band 
parameters may compete with the topological state. 16,18 
Moreover, local correlations related to the frequency de¬ 
pendent part of the self-energy can drive a topological 
phase transition. 19 These issues have been stimulated by 
research on iridium-based materials XilrO-$ (X = Na 
or Li) in which both spin-orbit coupling and electron 
correlations are strong. 20 Another new class of systems 


in which strong interactions could play important role 
are cold atoms in optical lattices. 21 Recent experimental 
progress in this field allows to test Hamiltonians with dif¬ 
ferent types of hoppings on lattices, synthetic spin-orbit 
couplings and synthetic gauge fields, as well as the inter¬ 
actions, in terms of their topological properties. 22,23 

Despite of intensive research, the influence of the 
strong correlations on the TIs at finite temperatures ap¬ 
peared so far only in few studies. 24-2 ' A major problem 
was how to extend a idea of the topological insulators 
to nonzero temperatures. Main progress in this field was 
made in Ref. [28] were a concept of topological phases of 
density matrices was introduced. However, it is still not 
know how the correlations affect the behavior of the TIs 
and how to compute topological invariants of interacting 
systems at nonzero temperatures. 

In this work we investigate the phase diagram of the 
correlated TIs at finite temperatures. Using the concept 
of topological phases for density matrices, we introduce 
methods to compute topological invariants for noninter¬ 
acting and interacting systems at nonzero temperatures. 
Then, we consider the Kane-Mele (KM) model without 
the Rashba spin-orbit coupling term. 1 To include corre¬ 
lations we extended the KM model by adding an interac¬ 
tion term as in the Falicov-Kimball (FK) model. It de¬ 
scribes a short-range local Coulomb interaction between 
itinerant electrons and localized spinless fermions. 29 We 
employ the Hartree mean-field and Hubbard I approx¬ 
imations to study the phase transitions and obtain the 
phase diagrams. This approach allows us to obtain semi- 
analytical results for topological phase transition and ex¬ 
amine the influence of the electronic correlations at fi¬ 
nite temperatures. Optical lattices offer the possibility 
of physical realization of this model. 21 A similar model 
has been investigated in Ref. [30] within the dynami¬ 
cal mean-field approximation (DMFT). 31,32 However, the 
long-range ordered phases have not been discussed in gen¬ 
eral cases. 30 

The paper is structured as follows. In Sec. II. we 
present the extension of topological invariants to finite 
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temperatures. In Sec. Ill we introduce the KM model 
with the FK type interaction. We present the condition 
for a change of a topological invariant in the system in 
Sec. IV and in Sec. IV B we describe the Hartree and 
Hubbard I approximations. In Sec. V we examine ther¬ 
modynamic phase diagrams with topological phases. In 
Sec. VI we present our conclusions. 


II. TOPOLOGICAL PHASES AT NONZERO 
TEMPERATURES 

A. Topological phases of the density matrices 

At a nonzero temperature the system is no longer in 
a ground state and there is a general problem how to 
define and interpret topological invariants. When the 
electron correlations are absent, the topological phases 
of insulators at zero temperature can be defined as ho- 
motopy equivalence classes of mappings k i —> H (k) , from 
Brillouin zone (BZ) torus T d to manifold of the gapped 
Bloch Hamiltonians (BHs). 8 ’ 33 For symmetry protected 
TIs, points of the manifold represent preserving the sym¬ 
metry gapped BHs. 

At a nonzero temperature a quantum system is de¬ 
scribed by a density matrix p of a mixed state. For 
translational invariant systems without electron correla¬ 
tions the density matrix takes the form p = Ilkesz p(k), 
where p(k) is called the Bloch density matrix (BDM). 
In order to extend the notion of topological invariants 
to nonzero temperatures and non-equilibrium systems, 
topological phases for the BDMs has been defined by us¬ 
ing an ensemble of nonorthonormalized pure states, c.f., 
Ref. [28], 

We express this definition in the equivalent but more 
convenient for our purposes form, which is analogical as 
for the BHs, i.e., the homotopy equivalence relation of 
the BDMs with a constrain that the spectral gap be¬ 
tween the valence band eigenvalues p v (k) of p(k) and 
the conduction band eigenvalues p c ( k) cannot be closed 
and protecting symmetries must be preserved during a 
homotopy transformation. Formally, 

Definition 1. let p( k) and q(k) be Bloch density matri¬ 
ces of size N of a d-dimensional system. Additionally, 
let corresponding to them density matrices be invariant 
under a symmetry group G. They are said to be topolog¬ 
ically equivalent, written as p( k) ~ q(k), if there exist a 
continuous map F : T d x [0,1] K > Dj^(N), such that: 

(i) V k {F( k, 0) = p( k) A F( k, 1) = 9 (k)}, 

(ii) Vt e [o,i] V k {pv( k,f) - p c (k,t) > 0}, 

(in) V seG V te[0jl] [F(t),U(g)} = 0, 

where F(t) is the density matrix corresponding to the 
BDM F(k,t). The eigenvalues p v ( c ){ k, t) are defined as 
n (n + l)th largest eigenvalue of the F(k, t), where the 


number n is determined by a physical parameter describ¬ 
ing occupation of states. In case of equilibrium systems 
it is the chemical potential p. It is a generalization of 
the definition for the BHs, where the valence and con¬ 
duction bands E v ( c \ (k) can be defined as n (n+ l)th low¬ 
est eigenvalue of H( k). The operator U(g) describes a 
symmetry operation g £ G. The space of the BDMs is 
labeled as Dn(N) = GL(C,N)/U(N). It is the manifold 
of strictly positive matrices of size N x N, because topo¬ 
logical properties are insensitive to scaling 28 and thereby 
the normalization condition of BDMs can be omitted. 

The Def. (1) applies to general BDMs and therefore de¬ 
scribes the topological phases for both thermal and non¬ 
equilibrium states. 34,35 The crucial role in the Def. (1) 
plays the spectral constrain ( ii ). From it results that 
to change the topological phase it is necessary to close a 
spectral gap or break a protecting symmetry. As a conse¬ 
quence of that, if two topologically nonequivalent systems 
belonging to the same symmetry class are connected, on 
the interface between them the transformation F( k. t) 
evolves through region of the space Djy(N) where the 
spectral gap disappears and thereby edge states occur. 
These states are insensitive to small perturbations pre¬ 
serving the protecting symmetry. Therefore, the spectral 
constrain (ii) establishes the bulk-boundary correspon¬ 
dence as in the case of BHs. 

Without the spectral constrain (ii) we can always de¬ 
form by means of a transformation F(k, t) all BDMs to 
the maximally mixed one and in that case there exists 
only the topologically trivial phase. The spectral con¬ 
strain is not present when all possible states all fully oc¬ 
cupied. This entails that sum of topological invariants of 
all bands must be zero. 


B. Topological invariants at finite temperatures 

We show that the Def. (1) provides for every finite tem¬ 
perature the same classification of the topological phases 
and formulas for the topological invariants as for BHs. 

1. Noninteracting systems 
In equilibrium state the BDM takes the form 

p -0H(k) 

P(k) = —(1) 

where Z M = Tr(e~^ H ), H is full Hamiltonian, /3 is the 
inverse of the temperature T and M is number of unit 
cells. Hence, for every T, Eq. (1) establishes one-to-one 
correspondence 

<p : H(k) p(k) (2) 

preserving the protecting symmetries and the spectral 
gap, i.e., p c ( k) = p v ( k) if and only if E c ( k) = E v ( k). 
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As a consequence of that, if two insulators belong to 
the same topological class of BHs, they must be in the 
same topological phase of BDMs and thereby p is iso¬ 
morphism of equivalence classes of topological insulators. 
Therefore, possible topological phases of finite tempera¬ 
ture BDMs have the same classification as in the periodic 
table of topological insulators and superconductors 6 and 
topological invariants are computed by substituting 

<^ _1 (p(k)) = -iln(Zp(k)) (3) 

into the place of H( k) in the standard formulas for the 
topological invariants. 6 ’' In the limit T —> oo this theo¬ 
rem is not valid because p is not injective function, i.e., 
does not preserve distinctness of elements. In that case, 
the system is always described by maximally mixed den¬ 
sity of states and thereby it is in the topologically trivial 
state. 


2. Interacting systems 

There is not known a general classification of topolog¬ 
ical insulators for interacting systems. Their topological 
invariants are constructed as an analytical continuation 
of the formulas for noninteracting systems. 8,36 For inter¬ 
acting systems the Hamiltonian cannot be expressed as 
a single particle operator. Therefore, in order to perform 
the analytical continuation of the noninteracting topo¬ 
logical invariants, they must be formulated in terms of a 
single particle Green’s functions G, 36 i.e., 

G _1 (w,k) =uI-H(k) (4) 

where ui is the frequency and I is unity operator. This 
generalization of topological invariants describes properly 
topological phases until the interacting system is contin¬ 
uously connected in the space of parameters of the model 
to some noninteracting system without closing the spec¬ 
tral gap. Hence, this description cannot be applied to 
fractional topological insulators. 

With additional assumption that T is treated on the 
same footing as parameters of the model, the above con¬ 
struction can be generalized to the density matrices. This 
entails that the interacting system can be adiabatically 
connected to a noninteracting one with a different tem¬ 
perature, what illustrates the line AB in the Fig. 1. 

As it was discussed in the Sec. IIB 1 all formulas for 
the noninteracting topological invariants are identical for 
BHs and BDMs due to isomorphism from Eq. (2). There¬ 
fore, by substituting Eq. (3) into Eq. (4) in place of 
H( k) topological invariants expressed by G(uj, k) can be 
calculated for the BDMs. The analytical continuation to 
interacting systems does not break this correspondence 
and thereby the same statements hold for the topologi¬ 
cal invariants of the interacting systems. In particular, 
also the concept of topological Hamiltonian 44,45 remains 
unchanged for the density matrices. 



x. 


FIG. 1. Exemplary topological phase diagram in space of 
interaction U, temperature T and parameter A. The red sur¬ 
face is the phase boundary separating the normal insulator 
(NI) and topological insulator (TI). The curve AB connects 
adiabatically two points. 


C. Interpretation of topological invariants at finite 
temperatures 

In recently published papers Ref. [37] and [38] differ¬ 
ent extension of the concept of topological invariants to 
finite temperatures is proposed. Those topological in¬ 
variants do not match with these derived here from the 
Def. (1). In contrast to our proposal, those topological in¬ 
variants depend explicitly on a temperature for noninter¬ 
acting systems. Additionally, it is shown in Ref. [28] that 
they are not unique. Another advantage of the Def. (1) is 
that due to the bulk-boundary correspondence the spec¬ 
tral constrain ( ii ) links our topological invariants to the 
physical observable, which is the number of edge states. 

Finally we note that at nonzero temperatures, the Hall 
conductivity is not quantized due to a thermal activation 
of states from a conduction band. Therefore, the topolog¬ 
ical invariants cannot be measured directly as response 
functions. However, at sufficiently low temperatures the 
contribution from the edge states should dominate. On 
the other hand, the edge states are experimentally acces¬ 
sible by angle-resolved photo-emission spectroscopy even 
at the room temperature. 10 


III. MODEL 

As discussed in the Sec. II, topological invariants for 
noninteracting systems at finite and zero temperature are 
the same, when parameters of the model are not tem¬ 
perature depended. Situation is more complex in case 
of correlated systems because the temperature change of 
the order parameter affects the self-energy of the elec¬ 
trons and hence the topological properties of the sys- 
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tem. To investigate this effect we use the KM model. 1 It 
is represented by the Hamiltonian which eigenstates ex¬ 
hibit the quantum spin Hall effect. 39 Formally, it is con¬ 
structed from two copies of the Haldane model 40 with 
opposite signs of the Peierls phase for the electrons with 
the up and down spins, respectively, moving on a honey¬ 
comb lattice. We extend this model by adding localized 
spinless fermions which locally interact with the itinerant 
electrons as in the FK model. 29 The full Hamiltonian is 


H = -t i cj a c ja -t 2 e#y 4 x'C' c . 7 V 

(i,j)c r <i,j»crcr' 




( 5 ) 


where Ci„ (c\ a ) is the annihilation (creation) operator for 
the itinerant electron with the spin a =\ (J,), and fi 
(//) is the annihilation (creation) operator for the local¬ 
ized spinless fermion at the lattice site i. The parameter 
t\ > 0 is the nearest-neighbor hopping whereas t 2 > 0 is 
the amplitude of the next-nearest-neighbor hopping. We 
use t\ = 1 which sets the energy units. The Peierls phase 
is (f>ij = for the next-nearest-neighbor hopping in the 
clockwise and anticlockwise direction, respectively. The 
third Pauli matrix is r z a , and it describes the change of 
the sign in the Peierls phase for the electrons with op¬ 
posite spins. The strength of the local FK interaction 
is given by U . In this work we consider a half-filling, 
i.e. one particle per unit cell for both kinds of fermions. 
Then the chemical potentials for the itinerant (/z c ) and 
localized (/if) particles are /i c = p/ =0. The Hamilto¬ 
nian is invariant with respect to: the lattice inversion, 
the spin rotation, the time reversal symmetry, and the 
particle-hole transformation. Due to particle-hole sym¬ 
metry the system does not possess an indirect band gap 
for any parameters. 

The localized particles in the FK model are thermo¬ 
dynamically coupled to the itinerant electrons and their 
spatial distributions is determined by a minimum of 
the thermodynamic potential of the system. On a bi¬ 
partite lattice at half-filling the solution of the Kane- 
Mele Falicov-Kimball (KMFK) model Eq.(5) with t 2 = 0 
possesses a long-range order at low temperatures, 41 in 
which the localized fermions form a checkerboard pattern 
and the itinerant electrons form a charge density wave 
(CDW) . Then the inversion symmetry of the ground 
state is broken. The transition between the CDW phase 
and the homogeneous phase is continuous. The order 
parameter is d = (n B — tia), where n a is the average 
number of localized particles per unit cell in the sublat¬ 
tice a = A or B. Increasing the hopping amplitude t 2 can 
lead to a change of the ground state. 42 We estimate that 
the CDW phase is stable for t 2 < 0.5t for an arbitrary U 
(see Appendix A). Therefore, we only present results for 
this regime of the t 2 hopping parameter. 


IV. TOPOLOGICAL PHASE TRANSITION 
A. General conditions 


Topological phases are quantified by topological invari¬ 
ants. For the model (5) its topological properties are fully 
determined by the topological invariant equivalent to the 
Chern number C for a single spin subsystem. 43 In the 
noninteracting limit and t 2 > 0 the system is in a topo¬ 
logically nontrivial phase with two edge states for open 
boundary conditions. 1 

In order to study topology of the solution for the model 
(5) with finite U we introduce the one-particle Green’s 
function G a (u> , k; d), which depends implicitly on the or¬ 
der parameter d for the CDW phase. The Green’s func¬ 
tion obeys the Dyson equation 

G“ 4 (w, k; d) = (u + fj, c )I - H 0 cr (k) - S CT (w, k; d), (6) 


where 


£<r(w, k; d) 


£cr a(u, k; d) 
0 


° 

T, aB (u, k;d) J 


is the diagonal self-energy matrix of the interacting sys¬ 
tem and fJo<r(k) is a matrix representation of the non¬ 
interacting Hamiltonian. From the theorem discussed in 
the Sec. IIB results that, topological invariant at every fi¬ 
nite temperature for our interacting system is determined 
by using the standard method 6 applied to the topological 
Hamiltonian defined by the noninteracting Hamiltonian 
matrix and the self-energy at u = 0, 44,45 i.e. 

K op ( k; d) = H 0a ( k) + £ CT (u; = 0, k; d). (7) 


Introducing symmetric and antisymmetric combinations 
of the self-energies 

^ ^ ... E CTj4 (0, k; d) ± £ ct b(0, k; d) 

^aS/RK*-, d) ^ 1 

the topological Hamiltonian for the spin er particles reads 

K° p (k; d) = H 0a ( k) + £ ctS ( 0, k; d)I + £ ctR ( 0, k; d)r z . 

( 8 ) 

The last term breaks explicitly the inversion symmetry. 
In rest of the paper we use local approximations in which 
the self-energy does not depend on the wave vector k, i.e. 
£ cr S/R(k;d) = E cr g/ R (d) is only a function of the order 
parameter d. 

Occurrence of the topologically nontrivial phase is con¬ 
nected with the fact that we cannot choose a one global 
gauge for eigenvectors of the Hamiltonian matrix (8) that 
is continuous and single valued over whole BZ. 46 This 
yields the condition 

|Sfl(d)| > 3V3f2 (9) 


for the transition from the topological phase to the trivial 
phase (see Appendix B). We note that this is a general 
condition for existence of nontrivial topology for any local 
approximation. 
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B. Considered approximations 

1. Hartree approximation 

Within the Hartree approximation the self-energy (6) 
of the itinerant particles is given by 47 

= ( 10 > 


Cf. Fig. 2. Since the electronic correlations are included 
within the Hubbard I approximation we conclude that 
they reduce strongly the topologically nontrivial solution. 
In details, this is due to formation of additional band 
states inside the Hartree gap, as is shown in Fig. 3, be¬ 
cause of the frequency dependence of the self-energy (13). 
This results that the energy gap is closed for smaller U 
in the comparison to the Hartree approximation case and 
thereby the topological invariant is changed. 


Hence, 


5Wd) = ~ (11) 

and from Eq. (9), the condition for the change of the 
topological invariant takes the form 


U>Uf 


6y/3t 2 

~w 


( 12 ) 


The Hartree approximation corresponds to the first order 
perturbation expansion with respect to small U. 


2. Hubbard I approximation 



Within the Hubbard I approximation the self-energy 
(6) of the itinerant particles is given by 47 


£<r(w,k; d) 


U n A 


U 2 n.An B 

u+(/(j-nB) 


U ns + 


0 'l 

U 2 nnn a I 

u]+U(^~n A )J 


(13) 


Hence, 


^an(d) = 


Ud U{ 1 - d 2 ) 
~2 2d 


U_ 

2d' 


(14) 


There are two contributions to Eq. (14): (i) the Hartree 
term, which takes into account a formation of the long- 
range order; and (ii) the term that takes into account 
effects of local electron correlations. Substituting 
into the Eq. (9), the condition for the change of the 
topological invariant takes the form 


U > U™ = 6V3t 2 \d\. (15) 


FIG. 2. Topological phase diagram of KMFK model within 
the Hartree (blue line) and Hubbard I (red line) approxi¬ 
mations in the space of (d, U *) parameters, where: U* = 
t//6\/3t2, TI - topological insulator, NI - normal insulator. 
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The Hubbard I approximation becomes exact when we 
neglect hopping ti and t 2 amplitudes. 


FIG. 3. Density of states of KMFK model within the Hartree 
and Hubbard I approximations. We use: t 2 = 0.2, U = 1.5, 
and d = 0.8. 


C. Topological phase diagrams 

In both Eqs. (12) and (15) the critical interaction 
strength U c depends linearly on t 2 . For |d| = 1 the 
correlation term in (14) vanishes and £7™ = £7®. For 
|d| < 1 the critical interaction U^ increases with decreas¬ 
ing |d| whereas £7™ vanishes linearly with decreasing |d|, 


In a homogeneous phase with d = 0 the noninteracting 
model has a topological insulating phase. For any finite 
interaction U > 0 the Hubbard I solution gives a topo¬ 
logically trivial Mott insulator. It means that U = 0 is 
a singular point in this approximation. The transition 
from the topological insulator to the normal Mott insu¬ 
lator occurs without closing the bad gap and it is only 
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driven by local correlations. We find that in the Hub¬ 
bard I approximation —> 0) = oo for any U > 0. 

However, the analytical condition Eq. (15) remains well- 
defined if we take d —> 0 limit at the end. Then, we 
do not have to use the concept of the frequency domain 
winding number 48 to compute the topological invariant 
in the homogeneous phase. 


D. More advanced approximations 


In Fig. 2 we see that the Hartree and Hubbard I so¬ 
lutions diverge from each other when the homogeneous 
phase is approached. Mathematically, it is due to the 
fact that when the long-range order vanishes the correla¬ 
tion part of the Hubbard I self-energy becomes singular. 
Physically, it reflects the fact that these two approxima¬ 
tions are valid in different limits, as discussed earlier. To 
resolve that problem one has go beyond the perturbation 
theory because one cannot get the Mott insulator by any 
local non-renormalized perturbation expansion (NRPE) 
in U 49 Indeed, within NRPE we find the following self¬ 
energy for the homogeneous phase 


H<r(w, k;d = 0) 


U U/ 4 

2 uj — Aq(w) ’ 


where the hybridization function Ao(w) = c o — Gq 4 {u}) 
and Go(w) is a noninteracting Green function at a single 
lattice site and is given by the Hilbert transform of a 
noninteracting density of states G o(w) = f po(co)/(lo — 
e). From the particle-hole symmetry G o (0) = — i7rp(0), 
where p(0) = 0 in our system. Hence, A 0 (w — > 0) = 
—ioo and thereby at oj = 0 the self-energy reduces to the 
Hartree self-energy. 

The most advanced renormalized local approximation 
is the DMFT. 31,32 Our preliminary DMFT results 50 show 
that up to the critical interaction strength U c \ the hy¬ 
bridization function at u> = 0 is the same as within the 
Hartree approximation and the system is a topological 
insulator. Then, the metal-insulator transition occurs 
and the hybridization function has finite imaginary part. 
At the interaction strength C / C 2 > U c \ the gap is opened 
again and the systems becomes topologically trivial Mott 
insulator with the hybridization function at uj = 0 equal 
to zero as within the Hubbard I approximation. Hence, 
the Hartree solution can be interpreted as an upper bond 
and the Hubbard I as a lower bond on the value of t/ cl ( 2 ), 
respectively. This result is consistent with a work, 30 in 
which the DMFT equations where solved in the homoge¬ 
neous phase. 


V. THERMODYNAMIC PHASE DIAGRAM 

Stability of the CDW phase is determined by minimiza¬ 
tion of the free energy with respect to the order param¬ 
eter d at a fixed temperature T and given microscopic 


parameters U and t 2 . The free energy of the KMFK 
model is 51 

F{d,T) = —^~ f p{w\d) ln(l + e~^ u )du> 

1 ( 16 ) 
-I- -[«a ln(n^) + tib ln(ns)], 

where the first term describes a contribution from the 
itinerant electrons and is equivalent to the free energy of 
the noninteracting fermions but with a density of states 
p(ur, d) modified by the interaction. The second term de¬ 
scribes an entropic contribution from the localized par¬ 
ticles. The inverse of the temperature is denoted by 
fd = 1/T. The temperature dependence of the topologi¬ 
cal phase transitions is obtained by inserting the optimal 
value d(T) into corresponding Eqs. (12) and (15). 
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FIG. 4. Phase diagram of KMFK model for t 2 = 0.2 within 
the Hartree approximation. 

The thermodynamic and topological phases within the 
Hartree approximation are displayed in Fig. 4 for a se¬ 
lected value of t 2 = 0.2. There are three phases: a ho¬ 
mogeneous topological insulator (TI-Hom), a topological 
insulator with the CDW long-range order (TI-CDW), 
and a normal insulator with the CDW long-range or¬ 
der (NI-CDW). The continuous thermodynamic transi¬ 
tion between the TI-Hom and the TI-CDW is marked 
by the red dashed curve. The topological transition be¬ 
tween the TI-CDW and the NI-CDW, displayed by the 
blue solid curve, is associated with closing a gap and a 
semi-metallic behavior right at the transition line. In¬ 
terestingly, the coexistence of the CDW long-range order 
and the topological phase occurs in a finite range of U 
values for a given temperature. This also contrasts to the 
Haldane model with nearest-neighbor interactions, where 
the charge-ordered phase is always topologically trivial. 52 
Since at T = 0 the order parameter d = 1 for any U, we 
find from Eq. (12) that the TI-CDW - NI-CDW transition 
takes place at U ® = 6\/3t 2 , which agrees with the Fig. 4. 
For other t 2 hopping parameters the phase diagrams are 
qualitatively similar to that presented in Fig. 4. How- 
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ever, when t 2 decreases the coexistence TI-CDW regime 
shrinks and finally disappears at £2 = 0 because it is t 2 
that induces the topologically nontrivial phase. 



U 

FIG. 5. Phase diagram of the studied model for £2 = 0.2 
within Hubbard I approximation. 

The thermodynamic and topological phases within the 
Hubbard I approximation are displayed in Fig. 5 for the 
same t 2 as in the Hartree approximation above. There 
are again three but different phases: a homogeneous nor¬ 
mal insulator (NI-Hom), TI-CDW, and NI-CDW. The 
homogeneous topological insulator, seen in Fig. 4, is re¬ 
placed by the homogeneous normal insulator. The topo¬ 
logical phase occurs together with the CDW long-range 
order and forms a bounded area on the phase diagram 
with a non-monotonic behavior of the critical tempera¬ 
ture. Again this is in contrast with the Haldane model 
with nearest-neighbor interactions, where the charge- 
ordered phase is topologically trivial. 52 We also note that 
at T = 0 the phase diagrams Figs. 4 and 5 are identical. 
This is a general feature of the solutions of the KMFK 
model within local approximations to the self-energy. 

The topological transition line in Fig. 5 has two charac¬ 
teristic behaviors: (i) for small U the thermodynamic and 
topological lines are almost tangent to each other; and (ii) 
for U close to U™ = 6a/ 3^2 the topological line is almost 
vertical. The course of the topological line results from 
Eq. (15) which relates the critical interaction strength 
I/™ and the value of the order parameter d(T) for which 
topological transition occurs. For temperatures close to 
the thermodynamic transition temperature T c , the order 
parameter takes the form d(T) ss (X\jT — T C (U), where 
a is a constant. Substituting this function into Eq. (15) 
the topological transition temperature T top is 

T «“ r ‘< c, > + ioSir <17) 

Hence, the topological critical temperature tends to the 
thermodynamic one, when U tends to zero and its main 
correction is of the order 0(U 2 ). In the second case the 
topological line follows the course of the order parameter 


when its value approaches unity. At low temperatures 
the order parameter is almost a constant, as is seen in 
Fig. 6, and this behavior gives rise to the vertical slope 
of the topological transition line in Fig. 5. 



FIG. 6. Order parameter as a function of the temperature for 
selected t 2 and U within the Hubbard I approximation. 

The topological phase diagrams for different t 2 are dis¬ 
played in Fig. 7. When the hopping parameter increases, 
an area of the topological phase is extended. This origi¬ 
nates from the fact that the hopping t 2 triggers the topo¬ 
logical phases in the KMFK model (5). The height of the 
nearly vertical part of the topological line decreases be¬ 
cause the constant part of the function d(T) decreases 
with increasing t 2 , as is seen in Fig. 6. 



u 

FIG. 7. Topological phase diagrams of the KMFK model for 
selected t 2 within the Hubbard I approximation. 

Comparing Figs. 4 and 5 we see that the critical tem¬ 
perature for the thermodynamic transition is reduced by 
the correlation effects taken into account in the Hub¬ 
bard I approximation and becomes non-monotonic. From 
Fig. 8 we see that the temperature of the thermodynamic 
transition decreases with increasing t 2 for each U. 
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Appendix A: Stability of CDW phase 



U 

FIG. 8. Thermodynamic phase diagrams of the KMFK model 
for selected within the Hubbard I approximation. 


In the strong interaction limit the KMFK model is 
mapped onto the antiferromagnetic Ising model, where 
Si = rii — 1/2 is a spin variable at a site i and an 
exchange coupling constants are J\ = t\/AU for the 
nearest-neighbors and J 2 = f 2 /4t/ for the next-nearest- 
neighbors hopping, respectively. 53,54 On the honeycomb 
lattice the ground state is a simple antiferromagnet for 
J2/J < 0.25 and a super-antiferromagnet otherwise. 55 
Such spin configurations correspond to the checkerboard 
and stripes pattern of the localized particles, respectively. 
Hence, we obtain the threshold ratio t^ltx < 0.5 of the 
CDW phase stability for arbitrary U. When it is not sat¬ 
isfied the CDW phase can still exist but only in a finite 
range of U. 

Appendix B: Change of the topological invariant 


VI. CONCLUSION 


In this paper we studied correlation effects on the 
phase diagram of the topological insulator at finite tem¬ 
peratures. We started with the definition of topologi¬ 
cal phases of density matrices and using it we extended 
a concept of topological invariants to nonzero tempera¬ 
tures. We introduced a model which allows existence of 
states with a long-range order, due to the interaction, 
and with nontrivial topology, due to a complex hopping 
amplitude. Using the concept of the topological Hamil¬ 
tonian and the dynamical local approximation 31,32 to 
many-body Hamiltonian we found analytically the con¬ 
dition for the topological phase transition. Within the 
Hartree and the Hubbard I approximations the complete 
phase diagrams with normal, topological, and long-range 
order phases were determined. The comparison of ther¬ 
modynamic phase diagrams allows us to conclude that 
both the temperature and the interaction play a crucial 
role for the existence of nontrivial topological states and 
they should be included in realistic calculations. In par¬ 
ticular, within the Hubbard I approximations, where the 
correlation effects are taken into account, the topological 
state appears only inside the CDW phase and is bounded 
to a finite area on the phase diagram. Moreover, the ho¬ 
mogeneous phase is a topologically trivial Mott insulator. 
The Hubbard I approximation seems to be the simplest 
approximation that can describe complete destruction of 
the topological phase only due to local correlations. 
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The topological invariant of the KMFK model (5) is 
equivalent to the Chern number for the spin up copy of 
the Haldane model. 43 Chern number is defined by 6 

C = —— / H(k)dk (Bl) 

2 71 " Jbz 

where the surface integral is over the Brillouin zone (BZ), 
and fl(k) = Vk x A(k) is a Berry curvature. The vec¬ 
tor potential A(k) = —i(u(k)|Vk|u(k)) is called a Berry 
connection and |u(k)) are eigenvectors of the lower block 
band of the Hamiltonian matrix (8). From Stokes the¬ 
orem we see that Chern number is nonzero when one 
global gauge of vector potential the A(k) does not exist. 
Therefore, this observation yields the practical test of the 
existence of a nontrivial topology. 

In the KMFK model, the topological Hamiltonian ma¬ 
trix is given by two dimensional matrix which can be 
written in the form 


H top ( k) = H x {k)r x + H y {k)r y + H z { k)r z , (B2) 


where the Pauli matrices are marked by t 1 with i = x,y, z 
and the decomposition coefficients H z ( k) are functions of 
the momentum. The spin index cr is omitted. For the 
two dimensional Hamiltonian matrix the eigenvector of 
the lower band, i.e. the eigenvalue for a given k, is 


*/(k) 


JL ( HZ W ~ I#001 \ 

N 1 \H x (k)+iH v (k) J 


(B3) 


where |if(k)| = y/H x (k ) 2 + FD'(k) 2 + H z ( k) 2 and N 1 is 
a normalization constant. This wave function is singular 
when H x (k) = H v (k) = 0 and H z ( k) > 0. We can gauge 
out the singularity by performing a gauge transformation 
rt J (k) -A u : ( k)e l< ^ k * obtaining the eigenvector 


u H ( k) 


i/( k)e^ (k) = 


1 

iv 77 


-H x (k) + my(k) \ 
H z (k) + \H(k)\ J' 
(B4) 










9 


It is now singular when H x (k) = H y { k) = 
0 and H z { k) < 0. The vector d{ k) = 

(H x {k),Hy(k),H z (k))/\H(k)\ sets a map from the BZ 
into the two dimensional unit sphere. If d(BZ) covers 
both poles of the sphere there is no a single gauge in the 
whole BZ to avoid the singularities. 

For the KMFK model H x = H v = 0 only at the Dirac 
points K = (47r/3\/3a, 0) and K = — K, where a is a 


lattice constant. At these points H z ( K) = Sfl(d)—3v / 3t 2 
and H z (K ) = T,n(d) + 3^/3t 2 - Hence, the vector d(k) 
contains only one pole for k gBZ when the condition 

|S fl (d)| > |3v / 3t 2 | (B5) 

is fulfilled and therefore the system is in the topologically 
trivial phase. 
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